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Рассматриваются возможности приложения результатов математического моделирования магистральных 
трубопроводов к задачам предварительного расчета технологических и конструктивных параметров участков 
транспортировки природного газа. Предлагается подход сквозного расчета дискретизированных значений 
переменных состояния при построении блочно-модульных схем моделей при их программной реализации, 
что позволит решать задачи повышения эффективности и параметрической оптимизации технологических 
свойств и технологической управляемости таких объектов при построении и эксплуатации АСУ ТП. 
Ключевые слова: участок магистрального трубопровода, математическая модель, установившийся режим, 
технологические и конструктивные параметры, аппроксимация, дискретизация. 


Введение. Для эффективного управления магистральными потоками природного газа, а также 
диагностики и прогнозирования состояния транспортной системы особое значение имеет возмож- 
ность математического и имитационного моделирования объекта. Это, прежде всего, касается 
задачи моделирования динамических процессов транспортировки газа в реальном времени [1]. 
Но математические модели пространственно распределенных объектов управления очень сложны 
по своей внутренней структуре взаимовлияния переменных состояния. Поэтому при разработке 
численных алгоритмов решения задач газовой динамики возникает проблема выбора внутренне 
устойчивой аппроксимационной динамической имитационной математической модели таких объ- 
ектов [2, 3]. Это касается способа аппроксимации переменных состояния и параметров сплошной 
среды, описываемой дифференциальными уравнениями в частных производных, их дискретными 
аналогами, т. е. разностными уравнениями [4-6]. Пространственная дискретизация приводит к 
аппроксимации производных в системе дифференциальных уравнений, являющихся математиче- 
ской моделью распределенного объекта. Далее вводится и дискретизация времени как основа 
численного решения полученных дифференциальных уравнений в полных производных. Такое 
преобразование связано с вводом в задачу параметров разностной сетки — пространственно- 
временных шагов. 

При таком способе решения задачи моделирования возможны два подхода [4, 7, 8]. Один 
ориентирован на определение закона дискретизации и аппроксимацию временных и пространст- 
венных производных непосредственно в дифференциальных уравнениях в частных производных 
[4, 7]. Другой подход основан на построении некоторого дискретного физического аналога, при- 
мером которого является так называемая ячеечная модель [9]. При этом система дифференци- 
альных уравнений в частных производных заменяется системой дифференциальных уравнений в 
полных производных, а размерность модели возрастает пропорционально количеству ячеек ана- 
лога. Тогда становится возможным использование хорошо развитой библиотеки алгоритмов чис- 
ленного решения дифференциальных уравнений в полных производных. 

В первом подходе система дифференциальных уравнений в частных производных заменя- 
ется системой алгебраических уравнений большой размерности относительно значений во всех 
узлах получаемой сетки пространственной дискретизации. Эти системы могут быть решены либо 
явно, либо итерационными численными алгоритмами (неявно). Такой подход позволяет получить 
пространственно-временную аппроксимацию с полным обоснованием свойств полученной систе- 
мы — сходимость, устойчивость, аппроксимация. Второй подход является вычислительно-экспери- 


470 


Вестник ДГТУ. 2011. Т. 11, № 4(55) 








ментальным, так как основан на сопоставлении свойств математической модели, аппроксимиро- 
ванной по пространственной координате и свойств численного метода решения дифференциаль- 
ных уравнений, т. е. аппроксимационного решения уравнений модели по времени. 

Для решения задач моделирования динамических процессов наиболее приемлем второй 
подход, но его реализация связана с множеством проблем, главными из которых являются устой- 
чивость решения и универсальность полученной модели [1, 4, 9, 10]. При этом если проблема ус- 
тойчивости модели нашла частичное решение в работах Н.В. Кудинова и проф. Р.А. Нейдорфа, то 
проблеме универсальности и неограниченной параметрической расширяемости модели до сих пор 
не уделялось серьезного внимания. Практика применения имитационной модели к задачам 
управления в рамках АСУ ТП придает этой проблеме особое значение. 

Решение неординарной задачи построения внутренне устойчивой аппроксимационной ди- 
намической имитационной математической модели упомянутых выше объектов невозможно без 
анализа фундаментальных результатов их исследования, т. е. без аналитически выраженных ма- 
тематических моделей. Кроме того, такие модели необходимы в проектных, технологических, ди- 
агностических и других расчетах, связанных с выбором наиболее эффективных конструкции или 
технологического регламента [11, 12]. Для корректного решения такого рода задач аналитически 
построенную и, следовательно, структурно адекватную процессу математическую модель необхо- 
димо идентифицировать параметрически. 

Структура взаимозависимостей переменных состояния распределенного объекта такова, 

что его математическая модель очень чувствительна к большинству конструктивных и технологи- 
ческих параметров и имеет физический смысл только в определенных, очень узких, их диапазо- 
нах. Исследования по определению этих диапазонов (как технологических, так и конструктивных) 
позволяют уточнить (т. е. параметрически идентифицировать) ряд трудноопределимых регла- 
ментных технологических переменных, а также рассчитать диапазоны их возможных значений. 
Кроме того, целесообразно разработать математические инструменты решения задачи выбора 
робастного сочетания конструктивных параметров [9]. 
Постановка задачи. Для универсальности аппроксимирующей распределенный объект имита- 
ционной модели необходимо придать ей структуру такой динамической системы, в которой мак- 
симально облегчается структурная перестройка, сохраняющая основные черты распределенного 
оригинала, но содержатся параметрические настройки технологических свойств, а также дейст- 
вия внешних возмущающих и управляющих воздействий. В связи с этим, проблему целесообразно 
исследовать на конкретном и максимально простом примере объекта с одной пространственной 
координатой. 

К таким объектам относятся трубопроводы транспортировки природного газа, которые в 
реальных магистралях реализуются в виде отдельных участков, разделенных газоперекачиваю- 
щими агрегатами. Такие участки магистральных газопроводов (МГП) и являются объектами 
управления в АСУ ТП. 

Решения, обеспечивающие блочно-модульное построение модели и алгоритма, эффектив- 
но раскрываются лишь на схемном уровне. Основой таких решений может быть только 
централизованное и сквозное обеспечение изменения всех настроек модели. Это относится как к 
конструктивным размерам, так и к технологическим параметрам. Кроме того, для обеспечения 
стабильности работы модели, необходимо в каждую ячейку программно передавать значения 
важнейших для объекта переменных состояния, отвечающих текущим статическим профилям их 
распределения. 

В связи со сформулированным подходом в статье ставится и решается задача расчета и 
сквозной передачи значений плотности и массового расхода газа. Наряду с этим важно получить 
оценки предельных возможностей исследуемого участка магистрального газапровода: конструк- 
тивных при заданных значениях технологических переменных и, наоборот, технологических при 
заданных значениях конструктивных параметров. Для решения всего комплекса обозначенных во 
введении задач необходимо исследовать те аналитические возможности, которые предоставляют 
получаемые на стадии математического моделирования динамические и статические характери- 
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стики распределенного объекта. Прежде всего, необходимо построить и оценить те взаимозави- 
симости переменных состояния и параметров объекта, которые позволяет выделить математиче- 
ская модель его установившегося состояния. 
Исходная математическая модель задачи. В исследованиях используется изотермическая 
математическая модель, задаваемая системой дифференциальных уравнений в частных произ- 
водных следующего вида, полученного путем нелинейного преобразования базиса переменных 
состояния модели: 

др_ 065. 


ры, 1 
01 Ох (и 


2 
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В (1), (2) введены следующие обозначения: р=р(х,!), и(х,!) и &==(х,1)= 


=р(х,1) *(х,г) - плотность и скорость транспортируемого газа, а также переменная состояния 


потока газа, имеющая физический смысл массового расхода, соответственно. Обе переменные 
являются функциями независимых аргументов: пространственного х и временного #. 

Константы модели и начальные условия переменных состояния заданы значениями, взя- 
тыми для одного из реальных участков магистрального газопровода: А =8,3144 Дж/моль.К; 


Т, =303К; и=0,01604 моль; ру =р(0, 0) = 47,7521 кг/м3; 5, = 506,268м’; Ор =1,22м; 


^. = 0, 000661 м/с. 

В выбранном для исследования виде математическая модель получена и исследовалась в 
многочисленных работах [1, 13-15]. Она является основой для получения различных частных мо- 
делей, изучение которых позволяет получить закономерности исследуемых связей и свойств объ- 
екта моделирования. 

Исследуемые математические модели (ММ) частных задач. Модель статики 
участка МГП. Эффективным приемом проверки корректности результатов построения динамиче- 
ской ММ является исследование ее в статическом состоянии. Для частного случая ММ (1)-(2), 
описывающего законы функционирования объекта в установившемс режиме транспортировки га- 
за необходимо учесть, что в соответствие с условиями статики 


Ор _0 
бр 04 0 
0 01 
Тогда из (1), с учетом отсутствия зависимости от аргумента 1, следует, что 
ЧЕ 
—=0->5=&, = сопзё, 
х 


и уравнение (2), с учётом однонаправленности векторной переменной = = &(х,1) в статике, при- 
водится к форме 

ар_  4Анрё» 

о Вбр)” 

Ах Б(ре, —АТр ) 
в которой переменные р и х разделяются: 


рр = А. (3) 


2 

РНЕ 
В связи с поставленной задачей необходимо проинтегрировать уравнение (3) и получить 
зависимость лишь одной переменной состояния р от х, так как в выбранном нелинейном базисе 


вторая переменная состояния & в статике становится константой. Исследуемый режим позволяет 
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проверить соответствие параметров стационарных точек фактическим их значениям, а также вы- 
явить взаимосвязи важнейших конструктивных и технологических параметров. 


Интегрирование (3) дает уравнение связи р с &, и х, которое оказывается неявным от- 
носительно р. В нем явно выраженным оказывается аргументх: 


р | 2 ря ())+ в (4) 
4. 255 Ро 

Эта особенность создает некоторые неудобства при расчетах, но позволяет, тем не менее, 
решать любые задачи, связанные с установившимся распределением плотности, а также через 


связь переменных состояния с распределением скорости потока 

==(х,1) = р(х,ё)\(х,!). (5) 
Влияние технологического режима транспортировки на конструктивные параметры 
участка МГП. Сложная по математической взаимосвязи переменных функция (4) параметриче- 
ски ограничена. Ее обратная структура задания аргумента по заданному значению функции по- 
рождает возможность попыток некорректного использования, когда подставляются недопустимые 
значения зависимой переменной и вычисляются абсурдные, с физической точки зрения, значения 
аргумента. 

В связи с этим функцию (4) необходимо исследовать на предмет получения реальных ог- 
раничений на область определения и область значений. Поскольку моделируемые участки не мо- 
гут иметь бесконечную протяженность в связи с протекающими при движдении газа диссипатив- 
ными процессами, а увеличивающаяся с падением давления и плотности газа скорость не может 
превышать скорость звука в его среде, функцию (4) необходимо исследовать на экстремум в сло- 
жившейся форме зависимости. 

Необходимое условие наличия экстремума 


ах пр ВТ Р( )+ 
= и 
2 
ар 4^| нЕ Р(х) 
приводит к следующему выражению, задающему граничное значение плотности газа при задан- 
ных значениях конструктивных параметров и технологических переменных: 


2 
6) 
х}=.|—. (7) 
Р(х) = рт 
Для определения характера экстремума длины (хотя он и очевиден) необходимо исследо- 
вать вторую производную функции (4) 
Ир аи 
ар’ 4\| на р’(х) 
Выражение (8) при любых р и 5, а значит и в точке (7), отрицательно, так что длина 


трубопровода, как и ожидалось, ограничена максимумом в этой точке. Таким образом, физикой 
процесса транспортировки обусловливается параметрическое ограничение 











=0 (6) 





<0. (8) 


2 2 
реле — а (9) 


НЕ. КТРо 

где Г, — длина исследуемого участка. 
Этот результат получен в работе [12] и исследован [13]. Здесь он иллюстрируется вспомо- 
гательными графиками. На рис. 1 представлена зависимость предельно допустимой длины участ- 
ка (Ат) от начального значения плотности (рии!). Формула и график позволяют оценивать 


как требуемую начальную плотность газа, обеспечивающую его транспортировку по участку за- 
данной длины, так и допустимую длину участка по номинальной выходной плотности газа вы- 
бранного ГПА. 
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Рис. 1. Оценка влияния плотности газа на выходе ГПА 
* 
на предельную протяженность участка МГП 


Аналогично представлена зависимость предельно допустимой длины участка ( Хи ) от ус- 
тановившегося номинального значения массового расхода газа (=пи/). Формула и график позво- 
ляют оценивать как допустимую длину участка по требуемой его производительности, так и пре- 


дельную производительность участка заданной длины (рис. 2). 


Характерно, что если в качестве технологического параметра брать установившуюся ли- 


нейную скорость потока газа м ( а] ‚ то предельно допустимая длина участка от начальной плот- 


ности не зависит. Это показано на рис. 3, где представлена зависимость предельной длины участ- 
ка ( Хэи ) от установившегося номинального значения скорости потока ( ипи/ ). Формула и график 
позволяют оценивать допустимую длину участка независимо от параметров газоперекачивающего 


агрегата, если известно требуемое значение начальной линейной скорости потока. 


* 
Здесь и далее формулы приводятся в записи и обозначениях, характерных для математического пакета Ма#САВ. 
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зх. 2 Хиц(10 ,47.75206) = 6661 х 108 
Хиц(18 47 75206) = 2056 х 108 
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Рис. 2. Оценка влияния требуемого значения массового расхода газа 
на предельную протяженность участка МГП 
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Рис. 3. Оценка влияния требуемого значения скорости потока 
на предельную протяженность участка МГП 
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Влияние конструктивных параметров участка МГП на его технологические свойства. 
При исследовании влияния конструктивных параметров на технологические удобнее строить об- 


ратные зависимости. Так, построена зависимость удельной производительности участка &(х,ро) 


от его длины х при различных начальных значениях плотности (р0) газа на выходе из газопе- 
рекачивающего агрегата (ГПА) (рис. 4). 


2 
Ес = т и Ес = 1.63361 х 10* 
2 у В(рО - ВС рб 


=[1 2.10747 75206] = 506.268 





9*10* 1х10° 1.1х10° 1.2х10° 1.3х10° 1.4х10° 1.5х10° 


х 
Рис. 4. Оценка влияния длины участка МГП на его удельную производительность 


Фактическая длина эксплуатируемого трубопровода Г, должна быть меньше предельной, 
а с учетом возможных изменений параметров участка МГП и проходящего через него потока, она 
должна быть меньше граничного значения с достаточным запасом. Так, например, рассчитанное 
для приведенных выше параметров участка МГП значение критической длины участка трубопро- 
вода составляет /,„, =258,3 км, т. е. реальный участок (120 км) почти вдвое короче предельного. 


Конструктивные параметры трубопровода рссчитываются на стадии проектирования. При 
эксплуатации существующей магистрали возникают вопросы, связанные с варьированием техно- 
логических параметров, в первую очередь, производительности. Полученные для статического 
режима зависимости позволяют проводить соответсвующие расчеты. Так, на рис. 5 приведена 
зависимость производительности участка МГП от начальной плотности газа рО при различных 


значениях длины участка. Расчет производится по формуле, полученной из (4) для полной длины 
участка Г: 





р ( 8 р? (Е))+ в?) (10) 


0 


=: р 
4). ео Ро 
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Рис. 5. Оценка влияния плотности газа на выходе ГПА 
на предельную производительность участка МГП 
при различных значениях его протяженности 


Выводы. Полученные результаты демонстрируют возможности, которые предоставляет инфор- 
мационная обработка статической модели участка МГП. При этом следует учесть, что рассматри- 
ваемый пример объекта и его математической модели упрощен за счет ряда допущений. Это, 
прежде всего, изотермический характер процесса транспортировки газа, равномерность распре- 
деления переменных состояния по сечению трубопровода, незначительность гравитационных со- 
ставляющих действующих сил, постоянство площади сечения трубопровода и др. Однако инфор- 
мационная насыщенность результатов такова, что позволяет делать обоснованные прикидочные 
расчеты для различных технологических задач и задач автоматизации. Использование адекват- 
ных моделей процессов транспортировки позволит реализовать все более тонкие и точные рас- 
четные технологии. 
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